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Abstract 

Two chains of Luttinger liquids coupled by both interchain hopping and 
interchain interaction are investigated. A degeneracy of two kinds of dominant 
states with in phase and out of phase ordering in the absence of the hopping 
is removed since the transverse fluctuation of charge density is completely 
gapful and that of spin density has two kinds of excitation with and without 
gap. The crucial effect of the interchain hopping on the electronic properties 
is studied by calculating susceptibilities. 
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I. INTRODUCTION 



It is well known that the ground state of one- dimensional interacting electron systems 
is non- Fermi liquid. The state, which is called as "Luttinger liquid", is characterized by a 
separation of spin and charge degrees of freedom, and anomalous exponents of correlation 
functions which depend on the interaction [jTj-Q. The parallel chains of the interacting elec- 
trons coupled by the interchain hopping are basic models of quasi-one-dimensional electron 
systems. Theoretical understanding of these systems is fundamental for studying electronic 
properties of organic conductors [Q, high temperature superconductors and interacting 
electron systems applied to strong magnetic field (Fractional Quantum Hall Effects) 0. In 
addition, the problem is important as a first step toward studying two-dimensional interact- 
ing electron systems. Consequently, there has been a growing interest in the system of the 
coupled chains, in particular two chains. Several work have been devoted to investigating 
the two chains coupled by the interchain hopping for the case of spinless Fermion ||^^12| and 
for the case including spin degree of freedom |[T3|-p2|. In addition, Anderson localization in 
such a system has been studied based on the above investigation 



However, interchain interaction also plays a role of the interchain coupling p5|j26[l . There- 



fore, by extending the work by Finkel'stein and Larkin |]T5| and that by Schulz |^, we study 
the two chains of Luttinger liquids in the presence of both the interchain hopping and the 
interchain interaction. To our knowledge, much is not known about such a problem in the 
case including spin degree of freedom. We clarify the electronic properties originated from 
the interchain hopping by calculating excitation spectrum, phase diagram and susceptibili- 
ties. The gap appears in the excitation spectrum of the transverse fluctuations. This leads 
to splitting of the degenerated states in the absence of the hopping, i.e., "in phase" ordering 
states and "out of phase" ordering states between the chains. From the calculation of charge 
and spin susceptibilities with and qy being the longitudinal and transverse wavenumber, 
it is shown that both susceptibilities with qy = ^ are the same as those in the absence of the 
hopping, and those with qy = n show remarkable dependence on q^. 
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The plan of the paper is as follows. In II, the Hamiltonian of the model is given and is 
expressed by use of the phase variables. In III, we study the excitation spectrum and the 
phase diagram at T = with T being temperature. The charge and the spin susceptibilities 
are also calculated. IV is devoted to discussion on the present results. 



II. MODEL AND PHASE REPRESENTATION 

A. Model Hamiltonian 

We investigate the system where two chains of Luttinger liquids are coupled by both the 
interchain hopping and the interchain interaction. The Hamiltonian is given by 

H = H]^ + Hint + Ti-'int ! (2-1) 

where 

= E efcp {4,p,a,i«fc,p,-,i + (1 ^ 2)} - t 5] {ai,p,^,iak,p,.,2 + (1^2)}, (2.2) 

k,p,a k,p,cr 

Hmt = ^ {(^ki,p,a,l(^k2-p,<7',l(^k2+q~p,a',iaki-q,p,a,l + (l ^ 2)| , (2.3) 

p,tr,a' k\,k2,q 

^(nt = ^^^E E {4i,p,<x,i«L-p,<x',2afc2+g,-P,-',2afci-g,p,a,i + (1^2)}. (2.4) 

p,(r,cr' k\,k2,q 

Equations (|2.2| ), ( p.3|) and ( p.4|) express the kinetic energy, the intrachain interaction and 
the interchain interaction, respectively. Quantities k and ekp{= vp{pk — kp)) denote the 
momentum and the kinetic energy of a Fermion where vp-, p = +{—) and kp are the Fermi 
velocity, the right-going (left-going) state of a Fermion and the Fermi momentum. The 
operator alp^^ expresses creation of the Fermion with k, p, a and i where cr = +(— ) and 
i{= 1, 2) are the spin I (|) state and the index of the chains. The interchain hopping and 
the length of the chain are defined by t and L, respectively. The normalized quantities, g2 
{g'2) denotes the matrix element of the interaction for the intrachain (interchain) forward 
scattering between particles moving oppositely. Note that the conventional definition of the 
elements are given hj g ^ g/{2-nvF) [0]- 



B. Phase Representation 



We represent the above Hamiltonian, Eqs.( p.2|) ~ ( p.4|) and the field operators of 



Fermions in terms of phase variables based on bosonization method. The separation of 
the Fermi wavenumber due to the hopping is taken into account by use of the unitary 
transformation, Cfc,p,o-,^j = {— f^o.k,p,cT,i + afc,p,o-,2)/v^ = ±). Then the kinetic term, Tit is 
rewritten as 

^k= 51 ^Fipk - kFf,)clp,a,^,Ck,p,a,^„ (2.5) 

where kp^ = kp — fit/vp. The interaction terms are rewritten as follows, 

T-Cint + ^I'nt = T^int + ^fnt + ^fnt' (2-6) 

'Hint = ^{92 + 92) J2 J2pp,^A(l)p-pyyi-<l)^ (2-7) 

p,<j,a' ,fj,,fj,' <? 

-^^nt = ^(92 - 92) E / H.,M^lp,.'>-P,<.',-M^P,<.,-4 ' (2.8) 
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p,o-,a-' ,fM 



n,3 _ T^^F 



2 {92-92) E /dx{^J^^^^^lp_^,_^z^_p,^/,/,z^p,<,,_^}, (2.9) 

P,a,CF' ,fJL 

where i}p^„^^ = {l^/L)J2kG'''''ck,p,a,^l■ The operator pp,a^^{q){= J2k cl+q,p,a,i,Ck,p,a,i,) expresses 
the density fluctuation around the new Fermi point kp^^ and satisfies the commutation 
relations, [pp,a,/.(-g), Pp',a',M'(5')] = 5 5 pp' 5 qq' 5 ^^^^ pqL / {2tx). The Hamiltonian Hl^^ denotes 
the interaction described in terms of the quadratic form of the density fluctuation. For the 
processes of Eqs.( p.8| ) and ( |2.9|) , the band index is not conserved in Tif^^^. and the index is 



exchanged in Tif^t, respectively. 



Here we deflne phase variables 6±{x), 0±(a;), 6±{x) and 0±(x) as p6|-pS 



e±{x) = 


F e(- 




(2.10) 


0±(x) = 




--\^\/'^^^^^j:a{p^,^,,{-q)±p^,^,,{-q)), 

(T,p. 


(2.11) 


~e±{x) = 


y e(- 


-"|^i/'+^''^^E/^(p+,^,.(-?)±p-,<^,m(-?)), 


(2.12) 


0±(x) = 


y e(- 


-1^1/2+^''^) a/z (p+,,,(-g) ± p_,,,(-g)) , 

4 


(2.13) 



where a ^ is a cutoff of the large momentum corresponding to the band width vpoi ^. 
Equations ( p.lO| ) ~ ( |2.13| ) satisfy commutation relations given by [6^{x),6^{x')] = 



[0+(a;),0_(a;')] = [e~+(.T), = = \n{l + [{x - x')a-^} - 

In {1 — i(a; — x')q;~^} ~ i7rsgn(x — x') and zero for the others. The variables, 6± and 
(f)± express the fluctuation of the total charge density and that of the total spin den- 
sity, respectively, while 6± and 0± express the transverse fluctuation of the charge den- 
sity and that of the spin density, respectively. Actually these properties are understood 

from the facts that dj+ = {T!-/V2)T,p,a,ii'l,a,ii^p,a,i, d^(j)+ = {TT/V^)T,p,a,iCri^l,a,ii^p,a,i, 

dj+ = (-7r/V2) Ep,a{^p,a,iV'p,a,2 + h.c.} and d^(j)+ = {-n/V2)J2p,a(^{'4'l,a,i^p,a,2 + h.c.} 
where ipp^a,! = i^/ VL)J2k^^''^'^k,p,a,i- In terms of the phase variables, the field operator of 



Fermions, ipp^a,ii{x), is expressed as p9| , |30 



= V^p,a,/.(a;)exp(i7rHp,^,^), (2.14) 



with 



©p,^,/. = ^{p^+ + ^- + ^'(p^+ + + + + ^f^ip^+ + ^-)}- (2-15) 

In Eq.( p.l4|) , the phase factor, vrSp o-.^t) is added so that the Fermion operators with different 
indices satisfy the anticommutation relation The factor Hp cr_^ is given by Hi = and 
Hj = I]}=i {i = 2 8) with iVj being the number operator of the Fermions with indices 
i where the index (p, cr, /i) corresponds to (+,+,+) = 1, — , -|-) = 2, (+,+,—) = 3, 
(+, -, -) = 4, (-, +, +) = 5, (-, -, +) = 6, (-, +,-) = ? and (-, -, -) = 8, respectively. 
Note that the above choice of Sp ^.^^ is not unique. 

By substituting Eqs. (|2.14]) and ( p.l5| ) into Eqs.( p.8|) and ( |2.9|) , and by using the fact that 

T^k = {T!-VF/L)J2p,a,tMT,qPp,a,f,{q)Pp,aA~l) PWl , the Hamiltoniau H = Ht + Hr can be 
expressed as, 

= |j / dx 1^(9.^+)' + Vp{dJ.f'^ + dx{{d,<P+Y + , (2.16) 



VF{g2 - 92) J I^Qg _ cog(2gox - V20+)} {cos v^0_ + cos , (2.17) 



+ 



where Vp = vf^{1 + 2g2 + 2g'^){l - 2g2 - 2g'^), r/, = ^{1 - 2g2 - 2g'^)/{l + 2g2 + 2g'^), 
Ag^± = A^^± = 1 and go = 2t/vp. Here Tix expresses fluctuations of total charge and 
spin densities, whose excitations are given by Vp\k\ and wfI^I, respectively. On the other 
hand, Tij?, which expresses the transverse fluctuation, includes the complex nonlinear terms. 
The nonlinear terms including cos \/26- result from 7i?^^ and those with cos(2goa^ — \/26+) 
are due to TYf^^. In deriving Eg. ( |2.17D , we chose a Hilbert space with the even integers for 



respective numbers Ni + N^, N2 + Ni, N^ + Nt, Nq + N^, Ni + N^, Ni + Nq and A^i + A^2 where 
Ni is the eigenvalue of iVj. The negative sign of cos(2goa; — y/2dj^) is due to such a choice 
of the phase factor, TcEp^cr^p. It is noticed that the interchain interaction does not change 
the structure of the Hamiltonian, i.e., the parameters in the presence of (73 are obtained by 
rewriting as g2 ^ 92 + 92 in "^t and 92 ^ g2 — 92 in T^R- 
The Hamiltonian, Tin can be rewritten as 



'Hr = vf j dx \%l)\{-id^)ipi - ^|(-i(9^)V'2 + ^^(-iSxOV's - V'4(-i<9x)^4} 
+ 'KVF{g2 - g'2) j dx - -44, - ^1^36-^'^"" - V^J^4e''^°"} 

X - i^2^i + + ^1^2 } , (2.18) 

where ■j/'j (z = 1 ~ 4) are new Fermion fields defined by 

ij, = ^=e'73(^++^-)e'f (^1+^2) = ^p[e'm+^^)^ (2.19) 
V 27ra 

ij2 = ^Le-^^('^+-'^-)e-'f (^^+^^) = V^^e-'i(^^+^^), (2.20) 
v27ra 

^3 = ^ 77^ (g+ +^"- ) (.i ^ (N3 +Ni)+in(N^ +N2 ) ^ ^'gi|(7V3+7V4)+i7r(JVi+^2)^ (2.21) 

v27ra 

^ ^ (.^i77^(g+-g-)(.-ii(jV3+jV4)+i7r(j^l+j^2) ^ ^/g-if (7^3+^4)+i7r(iVi+^2)_ (2.22) 

\/27ra 



In Eqs.( p.l9| )^( p.22"D , the phase factors with the number operator are introduced again to 



satisfy the anticommutation relation. Equation ( p. 181) is derived by choosing the Hilbert 



space with both A^i + N2 and + A^4 being even integers. It is noted that the phase 
variables in Eqs.( p.ll^ ) ~ ( p.22| ) are expressed as 



^~±(^) = -E^e(-"i''i/2+^'^^)(P3(-g)±P4(-g)) , (2.23) 
0±(:^) = -E^e(-"H/2^''^Upi(-?)±P2(-g)) , (2.24) 

where pj{q) = Eki'jik + q)^j{k) with = ^Y.k^Aky^'" (j = 1 ~ 4). 



III. PROPERTIES AT LOW TEMPERATURES 



By using renormahzation group method, Finkel'stein and Larkin [|T3| showed that the 



nonhnear terms in Eq.( p.l7[) with g2 7^ (^2 tend to the strong couphng in the hmit of low 
energy and the terms without (with) the misfit parameter, 2go become relevant (irrelevant). 
Schulz insisted that the transverse charge excitation is completely gapful and that of spin 
excitation has two kinds of excitation with gap and gapless from the symmetry of the 
Hamiltonian. Here we explicitly show the results by use of the mean field approximation 
in which the terms including the misfit parameter are neglected. This method is expected 
to be effective in the limit of strong coupling and has an advantage of the straightforward 
calculation of the several quantities. It should be noted that the break of the balance between 
6*+ and 9^ due to the misfit parameter may lead to the renormahzation of fjp defined by 
{Aq^^I Aq^^Y^"^ . We neglect such an effect as zeroth-approximation since the system tends 
to strong coupling and the gap appears. On the other hand, r]„{= (A^ ../A^ +)^/^) remains 
unity due to the balance between and 0_. The present method is effective for the energy 
lower than the hopping, and then the large momentum cutoff a^^ must be read as t/vp. 

A. Excitation Spectrum 

By making use of the mean-field approximation, Ti^ with (72 7^ 92 is rewritten as 

^i? = ^MF + ^MF - AA'L/(27rt;^((?2 - g'2)) , (3-1) 



where 

^MF = vf J dx - 'il>l{-id^)i>2} 

A/" 

+ - y dx {i^ui - i^2^i + i^Ui + ^1^2} , (3.2) 

^MF = VF [dx {4i-id,)i>s - i'li-id.)i'4} + A' [dx - . (3.3) 



The quantities, A and A' in Eqs.( p.2|) and ( |3.3|) are gap parameters determined by the 
following self-consistent equations, 

y = 7^vp{g2 - g',) {i (^M) - i (^4^3)} , (3.4) 

A' = nvpig^ - g',) {i - i + + {^l^} • (3.5) 
The eigenvalues, of Eq.(p]^) and 1^34 of Eq. (p.3|) are calculated as 



uji2= { (3.6) 



o;34 = ±^fe = ±V(^FA:)2 + A'2 . (3.7) 

The excitations of Eqs. (|3.6| ), whose spectral weights are 1/2, are obtained by T^mf terms 
of Majorana Fermion as 

n'^F = Y / dx {C^i-i^^)C^ - C2{-id,)C2} + iA J dxC^C2 

+ ^ 1 dx {Coi-id^)Co - Csi-id^)Cs} , (3.8) 

where Cq = i/v^(e^"/Vi - e"'"/Vi), C^i = 1/V2(e'"/Vi + e"'"/Vi), C2 = l/V2(e-^"/V2 + 
gi7r/4^t-j g^j-^^ _ i/y^(e~^'^/^-?/'2 — e^'^^'^ipl) ■ Finkel'stein and Larkin [jl5l have already ob- 
tained the same form as the first line in Eq. (|3.8| ) by replacing cos ^/26- in Eq. (p.l7|) with a 
value of a fixed point. However they did not show the gapless excitation which contributes 
to low energy properties e.g., specific heat 



The gap equations of Eqs.( p.4|) and ( p.5| ) are calculated (see Appendix A) as 



A = -2A'(^2 - 92) log ^ ^ , (3.9) 



f _|_ /£2 _j_ ^2 

A' = -A{g2 - g',) log ^^^g^ , (3.10) 



where is cut-off energy of the order of t. From the the gap equations, it is found that 
both (—A, —A') and (A, A') are solutions, and that sgn(AA') = —1 for g2 — g2 > and 
sgn(AA') = 1 for 5(2 - < 0- The solutions of Eqs.(pl9|) and ( p.lOD in the case of A' > is 
shown in Fig.|l|. 

B. Possible States 

We examine phase diagram which shows the most divergent state. Since the dominant 
contribution in the low energy limit is given by Tif^^ in Eq. (|2.8| ) [|15|, we rewrite the term as 

-{92-92) J2 / dx{(5I^lp,a,-M'^P,'^',A'')(II^i<x',pV'-p,a,-/.) 

-(^/^>lp,<,,_^/^/'p,a',^')(II/^V'p,a',/.V'-p,<7,-M)} • (3.11) 



p,(T,0" 

[92 - 92) 



Since the states should be selected so as to gain the energy from Eq.( p.ll|) , the possible 
states in the case from g2 — g2 > are given by 

ga,a' ^ J2f^^p,a,,,'ip-p,a',t, = " {^p,a,l^-p,a' ,2 + (1 ^ 2)} 

._ev. e.^' )*W )* sm|-^+p^^, + ^^_|, (3.12) 



^1 



7ra 1 72 2^2 2^2 J 

(3.13) 

and those in the case oi g2 — g2 < ^ are given by 

5+'"' = ^V^p,a,/.V^-p,a',/. = V^p,<7,iV^-p,<tM + (1^2) 

~ — ev^ e2v^^ ^'^+e2v^^ ' cos < ^ + p H i^4>-)i (3-14) 

TTQ! I V2 2v2 2v2 I 



1 -i2pkpx =^9+ ^{cT+a')<t>+ -^{a-a')^- / ^- , ^ ~ ^' 1 , ^ + ^' 1 \ 

7ra 1 72 2^2 2^2 J 

(3.15) 

Here DWl'" ( DW!^''^ ) expresses density wave with interchain and out of phase ordering ( 
with intrachain and out of phase ordering), while S'^'^ ( 5*+°^ ) expresses superconductivity 
with interchain and in phase ordering ( with intrachain and in phase ordering). The most 
dominant state between S'^'^ and DW^'" for (72 > g'2 {S'^" and DWl'" for g2 < g'2 ) 
is determined by the total charge and spin fluctuations. By noting that the correlation 
functions for the total fluctuations are calculated as, 

1 , 1 
(3.16) 

e^^+(^')e^^"+"'^'^+^"^e5^^"-"'^'^-^"^e^'+^°^e^^"+"'^'^+^°)e5^^"-"^ 

(3.17) 

we obtain the phase diagram in Fig.Q. The exponent of the correlation function of CDW is 
the same as that of SDW, and the exponent of singlet superconductivity is the same as that 
of triplet superconductivity, respectively. This comes from the symmetry of the interaction 
for the spin degree of freedom, as is also seen in spin dependent Tomonaga model with the 
isotropic interaction [Q. Therefore the phase diagram shown in Fig.^is essentially the same 
as that of spinless Fermion studied previously ||2^. In the repulsive case of 92 + 92 > 0, 
the most dominant state is the density wave. However, the interchain hopping leads to the 
superconducting state being subdominant even in such a region. 





C. Susceptibilities 

In this subsection, we calculate the charge susceptibilities, Xpilx, %',^<-^n), and the spin 
susceptibilities, X(t{Qx, %] ii^n), in the case of qx <^ Ikp and = or vr which are defined as 
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Xu{q..qy]ii^n)= dr J d{x - x')e''^"^e-"^^^^-^'hu{x - x' , Qy-, t), (3.18) 
and u = p or a. In Eq.( p.l8|) , 

Xp{x - x', qy- r) = ^ (t. {p{x, 1; r) + e^-^vl^;, 2; r)} {p(x', 1; 0) + e^'y p{x' , 2; 0)}) , (3.19) 

Xa{x - x', Qy; r) = ^ (Tr {m{x, 1; r) + e^'^^mix, 2; r)} {m(x', 1; 0) + e'«^m(x', 2; 0)}) , (3.20) 

where p{x, i; r) = J2p,a i'l,a,i{x; T)tpp^„^i{x; r) and m(x, r) = J2p,a ^)i'p,aA^'^ ^) de- 

note operators of the charge and spin densities at the i-th chain, respectively. 

At first we consider the case of Qy = 0. From Eg. ( |2.16D , both Xp{x — x',0]t) and 
Xaix — x', 0; r) are calculated as 

Xp{x - x', 0; r) = ^ {T^+ix, T)d,'9+{x', 0)) 

irvp (3L ^ ujr? + {vpq^^Y 
Xa{x - x', 0; r) = ^ (T^9a;0+(x, r)93;/0+(x', 0)) 

TTt;^ /3L a;„2 + (^;^g,)2 • ^ • ^ 

Therefore (itu^ — cj), one obtains 

Rexp(g., 0; c) = , (3.23) 



2 [vpqx 



\2 



Rex.(g., 0; cu) = , ^ , , (3.24) 

which are familiar to Luttinger liquid |]1|. 

Next we examine the case of qy = vr, for which Eqs.( |3.l9| ) and ( |3.2(JD are expressed as 

Xp{x -x',7r;T) 

= lv^\ 51 ^^^a^-/.'(a;;r)^/'p^<,^;,'(x;r) J ( ^ ^/^J,,_^(x0^p,<x,-M(a;'))), (3-25) 

\ \ p',(7',fj.' / \P,0-,M / / 

Xaix - x',tc;t) 

= lvA ^Vj',.',-;,'(a;;^)^p',a',/.'(x;r) J [ ^ (^■iplaA^')i^P'-^-p{x')]) ■ 

\ \p',(7',n' / \P,0-,/i / / 

(3.26) 
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After some manipulations (Appendix B), the static susceptibilities, ReXp('?a;, tt; 0) and 
ReXailx, vr; 0) at T = are respectively obtained as 

Rexp(g.,vr;0) = - / dJ (l - i^i^L + 

_^ \ /-j^ ^fc ik+q^-qg ^ A \ I (3 27) 

Ek + ^fc+g, -go Ek ^fc+g, _go Ek ^fc+g, -go J ' 



Rex.(g.,vr,0) = — Py dfc(l + 



E'J E'k-ik + VF{q. + qo) 



+P I dkll-^^ ^ 



Qx+qo V Ey Ej^ + ^k- vpiqx + go) 
r-qx+qo f Ck\ 1 



+P dkll- 



Qx+qo 



E'k) E'l, - + vpi-q. + qo) 
E'^ + ^k-VF{-q, + qo)j 



^ < -^J:,4V^4i.4)|. (3-) 



where P denotes the principal value, = vpk and e± = vpi^qx + '?o)- In Fig-i| and 
Fig.^, we show the normalized quantities Xpiqx,'^]^) and Xpiqx,'^]^) which are defined by 
ReXp(5'i, tt; 0)/(2/7rt>j7) and ReXaiqx,'^,0)/{'^/T^vp), respectively. The cutoff energy, ^c, de- 
fined in Eqs. ( |3.9| ) and ( ^.101 ) is taken as 2t. Note that Eqs. ( ^.27] ) and ( |3.28D are valid in the 
case of \qcc ± qo\<qo. 



Equation (|3.27|) shows that the value of ReXp{qx, tt; 0) in the case of g2 — g2 < is larger 



than that in the case of — > because g2 — 02 < {> 0) leads to sgn(AA') = 1 (—1). 
Such dependence on g2 — §2 is also found in the absence of the hopping (see Eq. ( |C14|) ). 
On the other hand, ReXo-(5'x, tt; 0) is independent of sign of the relative interaction. This 
result seems to correspond to the fact that ReXai^qx^T^]^) in the absence of the hopping is 
independent of g2 — §2 (see Eq. (|C15| )). 

In Fig.|, ReXp(g^, tt; 0) takes a minimum around q^ = qo in the case of g2 — g2 > 
and has a small dependence on q^ and g2 — g'2, i-e., being nearly unity for g2 — g'2 < 0. 
On the other hand, ReXo-(5'x, tt; 0) in Fig.^ takes the minimum for both g2 — g'2 > ^ and 
92 ~ g'2 < ^- The characteristic dependence is due to the separation of Fermi wavenumber 
i.e., kp- — kF+ = qo and the gap in the transverse fluctuation, both of which result from the 
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interchain hopping. Note that the g^:- dependence of these ReXpi^x, tt; 0) and ReXaiQx, tt; 0) 
does not change quahtatively by the choice of t- 

IV. DISCUSSION 

In the present paper, we studied the low temperature properties of two chains coupled 
by the interchain hopping and the interchain interaction where the interactions of only the 
forward scattering between oppositely moving particles were taken into account as a simplest 
model of Luttinger liquid. 

There are four kinds of excitations originated from the fluctuations of total charge, total 
spin, transverse charge and transverse spin respectively. The total fluctuations which show 
the gapless excitation are the same as those in the absence of the interchain hopping (Ap- 
pendix C) . On the other hand, the transverse fluctuations of both charge density and spin 
density which are expressed by the complicated non-linear terms are crucial in the presence 
of the interchain hopping. By utilizing the mean field approximation, it was shown that the 
transverse fluctuation of the charge is completely gapful and that of the spin has the two 
kinds of excitations with and without gap. 

The most dominant states are obtained as DWl''^ for g2 > \g'2\, DWZ''^ for g'2 > \g2\, 
S^i^" for g2 < — \g2\ and S'^'^ for g'2 < —\g2\7 respectively where density wave (superconduc- 
tivity) belongs to out of phase (in phase) ordering between the chains, i.e., the transverse 
wavenumber being 7r(0). We note that in the case of quasi one-dimensional electron system 
with only the hopping of pairs the density wave with the transverse wave number (vr, vr) 
or superconductivity with (0, 0) has maximum critical temperature. It is worth while noting 
that the states of the superconductivity remain subdominant even for the repulsive interac- 
tion, i.e., 5'2 + fi'2 > 0- Such a result may be the important point toward understanding of the 
competition of superconductivity and SDW observed in quasi one-dimensional conductors, 
(TMTSF)2X [0. 

Possible states in the absence of the interchain hopping are obtained by calculating 
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correlation functions for S'[[''^ (= ijp,a,ii'-p,a' ,i or 'ipp,a,2i^-p,a' ,2) , S'^'' (= i'p,a,iip-p,a' ,2 or 
'^p,a,2'ip-p,a',i), DW^^" (= V'p^^^iV'-p,^',! or ?/'^_^^2V'-p,<t',2) and DW^X''^ (= ^p,a,i^-p,a',2 or 
'^p,o-,2^-p,o-',i) which express the order parameters of the intrachain superconductivity, the 
interchain superconductivity, the intrachain density wave and the interchain density wave, 
respectively (Appendix C). The phase diagram in the absence of the interchain hopping is 
shown in Fig.|^. By comparing the phase diagram in Fig.|^ and that in Fig.^, it is found that 
the energy gain due to the interchain hopping removes the degeneracy of "in phase" and 
"out of phase" ordering. 

The g2;-dependence of charge and spin susceptibilities were calculated for both qy = Q 
and Qy = IT. The susceptibilities in the case of g^^ = are the same as those in the absence 
of the hopping since the total dynamics is not affected by the hopping. On the other hand, 
the static susceptibilities, ReXp('?x, tt; 0) for g2 — g'2 > and ReXo-l^z, tt; 0) in the case of 
Qy = n show the minimum around Qx = Qo which is ascribed to the separation of the Fermi 
wavenumber and the excitation gaps of the transverse fluctuation in the presence of the 
interchain hopping. The fact that ReXpiQx, tt; 0) in the case of g2 — g'2 < is larger than that 
in the case of g2 — g2 > is found also in the absence of the interchain hopping. 

We treated the system with the interaction processes of the forward scattering between 
the oppositely moving particles as the simplest model of Luttinger liquid. However, the two 
chains coupled by the interchain hopping have been known to show the various properties 
depending on the parameters of the system. The repulsive backward scattering becomes 
relevant in the low energy limit and opens the gap in the excitations of the total spin, and 



thus modifies the electronic properties |2^. This fact is different from the strictly one- 
dimensional case, where the repulsive backward scattering is renormalized to zero. The 
two chains of Hubbard Model shows the richer phases depending on the magnitude of the 
intrachain interaction, the interchain hopping and the filling [^^^. In addition, it has 
been reported that the two chains coupled by both the Coulomb repulsion and the exchange 
interaction show the superconductivity |^3|. Therefore further investigations are needed to 
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identify the ground state of two coupled chains as a crossover from one dimension to higher 
dimension. 
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APPENDIX A: GREEN FUNCTIONS OF EQS.(3.2) AND (3.3) 



We calculate Green functions of Eg s . (13 .21 ) and (|3.2| ). By using the solutions of Majorana 
Fermions which are calculated from Eq. (|3.8|) , the Green functions corresponding to the 
transverse spin fluctuations are calculated as 

1 



- (Tr4'l{x,T)i)l{x',T' 
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where is the time ordering operator, = vpk, 1/ [3 = T and = (2n + l)7rT. On the 
other hand, by diagonalizing Eq. (|3.3|) , the Green functions expressing the transverse charge 
fluctuations are calculated as 

1 
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(A3) 
(A4) 
(A5) 
(A6) 
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From these Green functions, one obtains in the limit of T = 



(A7) 
(A8) 
(A9) 



^2^1 ) = -i (V^2V^i ) = -i {i'li'2 



{iA/(4vrt;^)} \og{{^, + + A2)/|A|} 



^4V^3 = {iA7(2vrt;^)}log{(ec+ VC' + A'2)/|A'|} . 



(AlO) 
(All) 



By substituting Eqs.( [Al(j| ) and ( |A11D into Eqs.(|3.4|) and (P75|), one obtains Eqs.(^.9D and 
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APPENDIX B: CALCULATION OF SUSCEPTIBILITIES, EQS. (3.27) AND (3.28) 



Neglecting the constant whose absolute value is unity, the quantity -Fp(o-) 
J2p,aA^)'^l,<7,f,i'p,'T-^i is expressed as, 



1 



2Tia 



+ expi-^(0++_ 



(Bl) 



which is derived by calculating (T^fJ, r)Fp(o-)(a;'; 0) ) with the precise treatment of the 



negative sign originated from the phase factors in terms of Np^^^^ and iVj (i = 1 ~ 4) . Equa- 
tion (pip shows that the susceptibilities with Qy = tt are expressed by only the transverse 
degree of freedoms. 

By use of Eqn.(^) ~ (^), the quantity (^Vi^pV) (x; r)Fp(a.)(x'; 0)^ is calculated as 

(T,F;(^)(a;;r)F,(.)(x';0)) 
= e-'«»(^-^')(T. { [Mx; r) - (+)e-/2^I(x; r)] ^3(0:; r) - [4{x; r) - (+)e-/2^2(x; r)] r)} 

X {V^l(x') [V^l(x') - (+)e-'-/2^i(x')] - ^4(x') [V2(x ) - (+)e-^-/2^|(x')] } ) 
+ (x ^ x', r ^ — r) 



e-'«° i/p(,)(x - x', r) + (x ^ x', r ^ -r) 



(B2) 
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and the factors, m^, Vk, u'^, and are defined by 



Uk = 1 
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Vk = 1 
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Wfc = 1 
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v'k = 1 




(B8) 



From Eqs.(p3D and (@), xpi.Qxy'^A^n) and vr, iti;„) are calculated as, 
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(BIO) 



where f{z) = l/(e^^ + 1) is a Fermi distribution function. Equations ( [B9D and ( |B1(J| ) in the 
hmit of T = lead to Eqs.{^2^ and (^M )- 



APPENDIX C: PROPERTIES IN THE ABSENCE OF THE HOPPING 

In this case, by using density operators, pp^cr,u{(l){= J2k(^k+q,p,cr,u(^k,p,cr,u) with u = + and 
V = — being the chain index i = \ and i = 2 respectively, Eqs.( p.2|) ~ (|2.4| ) in the absence 
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of the hopping can be written as, 

p,o-,^ q 

Hint = ^^^y^ ^ ^ {pp,aA(i)p-P,'^A-(i) + pp,<rA(i)p-P,'<rA-(i)} > (C2) 

p,cr,u q 

'X-'int = ^^^7^ H H {Pp,a,z.(g)P-p,a,-i.(-g) + Pp,a,^(g)p-p,-a,-i.(-g)} ■ (C3) 

p,cr,u q 

Then the Hamiltonian is rewritten in terms of the phase variables as 

+ ^ / dx {(9.0+)^ + (4.0-)' + a.0V)2 + (4.0'_)'} , (C4) 



where 
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,^0 v^9^ 
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^'''-)E^Kp+,'^.^(-9)±p-,^,^(-?)) ' 


(C8) 



and ^, = J(l - 2(72 - 2(7^)/(l + 2g2 + 2g'^), 7]' = J{1 - 2g2 + 2(7^)/(l + 2^72 - 2(7^), v 



VF^il + 2(72 + 2(70(1 - 2(72 - 2(70 and v'^ = Vp^Jil + 2(72 - 2(70(1 - 2(72 + 2(70. The quan- 
tities, 6'± and 0± defined here are the same as those of Eqs.( p.lO| ) and ( p.ll| ), while 9'^ and 



0'^ are different from 9± and 0± of Eqs. (|2.12| ) and (|2.13|) . Since the parts including 9± and 
0± in Eqs.( p^) is the same as Eq.( p.l6|) , the dynamics of total fluctuation does not depend 
on the interchain hopping. By making use of field operators, ipp^a^u expressed as 

+ 9^ + u{p9'^ + 9'_) + (t(p0+ + 0_) + (Tz/(p0; + 0'_)} 



e'' 



lljp,a,v = r- exp 

V27ra 



2^2 



(C9) 

correlation functions of the order parameters are calculated as. 
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SrHx)S^^'^iO)) r^^-j , (CIO) 

SrHx)Sr{0)) ^^-^^j , (Cll) 

DW[l''''^{x)DW[l'''\0)) ~ f , (C12) 

(5r^(x)^r(o)) ~ (^^j , (C13) 

where S"^'"' = ^pp,a,,yi^-p,a',u, S'^'"' = Vp,<7,i.V'-p,a' -i., DW^'"' = V^J and DW^''"' = 
ip^^^^^ip-py-u express the order parameters of the intrachain superconductivity, the inter- 
chain superconductivity, the intrachain density wave and the interchain density wave, re- 
spectively . In Eg. ([C^) , we neglect the phase factor given by Np^a,u, because the particle 
number of each branch is conserved in the absence of the hopping. The exponent of the 
correlation functions of "in phase" ordering is the same as that of the "out of phase" order- 
ing in addition to degeneracy in spin degree of freedom. The phase diagram obtained from 
Eqs. dCToD ~ (|CT3|) is shown in Fig.|. 

Susceptibilities are calculated as follows. From Eq. (|C4|) , it is found that ReXpiqx,^,^) 
and ReXa{<lx,0,u!) are the same as Eqs. (p.23|) and ( p. 241) . By noting that 



Xp{x - x', ^,r) = — (T^dJ'^ix, T)dx:e'^{x', 0) 
Xa{x - x', Tr,T) = ^ (Trdx(j)'^{x, r)9x'0+(x', 0) 

we obtain ReXp(5'x5 tt, cu) and ReXaiQx,'^,^^) written as 

Rexa(g^.,vr,u;) = r, (C15) 

TlVp [VFqxj — ^ 

which are the forms familiar to the Luttinger liquid. Note that RexpilxiT^i^) is suppressed 
(enhanced) in the case of (72 — > (< 0) while ReXcr(5'x, tt, cu) is independent of (72 — 92- 
Thus it turns out that the two chains coupled only by the interchain interaction remain as 
the Luttinger liquid. 
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FIGURES 

FIG. 1. Solutions of the self-consistent equations, A/^c and A'/^c as a function of g2 — 92, 
which are obtained from Eqs.(3.9) and (3.10). Note that (—A, —A') is also the solution as well as 
(A, A'). 

FIG. 2. Phase diagram in the presence of the interchain hopping. Here DWj {DW^' ) 
and 5^'^ {S'^^ ) express density wave with interchain (intrachain) ordering with out of phase and 
superconductivity with interchain (intrachain) ordering with in phase, where a and a' express spin 
indices. 

FIG. 3. Normalized charge susceptibility, xp{<lx,'7^,0){= Rexp{qx,''^',0) /{^/t^vf)) at T = as 
a function of vfQx/^c in the case of g2 — 92 = -0.2, 0.2, 0.3 and 0.4 where = 2t. 

FIG. 4. Normalized spin susceptibility, Xa{(lx,'^',Q){= R-eXo-C^x, ti"; 0)/(2/7rvF)) at T = as a 
function of vfQx/^c in the case of \g2 — 52! = 0-2) 0.3 and 0.4 where = 2i. 

FIG. 5. Phase diagram in the absence of the interchain hopping. Here DW^'"^ {DW^'"^ ) and 
S l_ (^ii' ) express density wave with interchain (intrachain) ordering and superconductivity with 
interchain (intrachain) ordering, where a and a' express spin indices. 
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